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Abstract. Consider a worldline of a pointlike particle parametrized by polynomial 
functions together with the light cone (“retardation”) equation of an inertially moving 
observer. Then a set of apparent copies, R- or C-particles, defined by the (real or 
complex conjugate) roots of the retardation equation will be detected by the observer. 
We prove that for any “polynomial” worldline the induced collective dynamics of R- 
C particles obeys a whole set of canonical conservation laws (for total momentum, 
angular momentum and the analogue of mechanical energy). Explicit formulas for 
the values of total angular momentum and the analogue of total rest energy (rest 
mass) are obtained; the latter is “self-quantized”, i.e. for any worldline takes only 
integer values. The dynamics is Lorentz invariant though different from the canonical 
relativistic mechanics. Asymptotically, at large values of the observer’s proper time, 
the R-C particles couple and then assemble into compact incoming/outgoing clusters. 
As a whole, the evolution resembles the process of (either elastic or inelastic) scattering 
of a beam of composite particles. Throughout the paper the consideration is purely 
algebraic, with no resort to differential equations of motion, field equations, etc. 


1. Introduction 

We continue to elaborate the concept of the so-called “one-electron-Universe” initiated 
long ago by J.A. Wheeler and R. Feynman BIZI- The question is about the construction 
of self-consistent dynamics of a set of identical particles that all belong to one and the 
same worldline. Specifically, in M we considered the class of pointlikc particles’ 
worldlines defined implicitly by a system of algebraic equations 

^t a (x,y, z,t) = 0, a = 1,2,3. (1) 

where {'!'„} are three arbitrary independent functions of three spatial coordinates and 
the timelike parameter t. Then, for any value of t, one has a finite set of roots 
{xk, Vk, Zk}, k = 1,2, ...,1V defining the positions of N pointlikc particles. The 
latter are considered identical and possess, in particular, equal masses and charges. 
In the course of time, the particles-roots move along the trajectory curve resulting from 
(JT|) (through elimination of t) which, nonetheless, may consist of a number of isolated 
branches. Anyway, the fact that all the particles belong to the same “worldline” imposes 
rigid restrictions on their dynamics since the roots are strongly correlated. 

These correlations become manifest when the generating functions {Ta} in ([T| are 
arbitrary polynomials , both in {x a } and t. Indeed, after two of three spatial coordinates 
are eliminated, the system reduces to a polynomial equation for one variable, with 
coefficients being some polynomials in t. For such a system, the above correlations 
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are explicitly represented by the well-known Vieta formulas. Moreover, for any non¬ 
degenerate (in the sense specified in [3]) form of generating polynomials these formulas 
give rise to a system of conservation laws reproducing standard laws of Newtonian 
mechanics, namely, the law of conservation of total momentum and (the analogue 
of) total mechanical energy. Total angular momentum is also conserved for any non¬ 
degenerate polynomial system. 

Our main goal in the paper is to make an attempt to put the scheme in 
correspondence with the requirements of Special Relativity. The realization is 
undertaken in Sect.2 and 3 where we consider a “unique worldline” and an inertially 
moving observer, together with corresponding retardation equation (equation of the 
observer’s light cone) marking a set of points on the worldline. Note that such points, 
in fact apparent copies of one and the same particle, have been considered previously in 
[5] Ej (for tachionic worldlines) and in [8) [9] (for the case of the complexified space-time 
background). 

In the paper we consider only the familiar class of worldlines defined in the canonical 
parametric way, i.e. by X /Jr = i^(r), /i = 0,1,2,3. Even in this simplest case 
one obtains, from the principal equation of the light cone of an “inertial” observer, 
an ensemble of identical particles with nontrivial correlated dynamics. Again, the 
correlations become manifest when the parameterizing functions (i 7 ),} are arbitrary 
polynomials in r. In this case, we immediately get the polynomial equations for each 
of the coordinates {A^} and the parameter r itself, with coefficients depending on the 
(proper) time of the observer T. Applying then the Vieta formulas, one again obtains, 
for an arbitrary “polynomial” worldline , a full system of conservation laws for the set of 
roots-particles! 

Note also that one should take into account both real and complex conjugate roots 
of the considered polynomial system of equations. The latter have equal real parts and, 
w.r.t. to the latter, can be visualized in the observable 3-space. On the other hand, 
both parts of complex conjugate roots enter into the Vieta formulas and, consequently, 
contribute to the conservation laws. Thus, the complex conjugate roots should be 
identified with the second, composite kind of particles (with twice greater mass); in [3, 3] 
they were called C-particles (to distinguish with R-particles corresponding to real roots). 

The Lorentz invariant [] structure of the equation of light cone allows one to 
regard the whole scheme as a relativistic generalization of the previously developed 
Newtonian-like algebraic mechanics. Then one should consider as physically meaningful 
only the SO( 3, l)-covariant characteristics of the dynamical system of R-C particles. 
This requirement results in a rather specific set of the conservation laws which do not 
fully reproduce those specific for the canonical relativistic mechanics. 

Moreover, velocities of R-particles are, generally, not limited. Indeed, as in the 
purely Newtonian case pun, at some instants of the observer’s time T he observes an 
event of merging of a pair of R-particles when some two real roots become multiple 


1 In fact, one can exploit the more general conformal symmetry of the light cone equation 
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and then transform themselves into a pair of complex conjugate ones. Physically, this 
situation models the process of annihilation of two R-particles (one of which should 
be considered as an antiparticle) accompanying by the formation of one composite C- 
particle. At the points of annihilation, singular points of the worldline, velocities of the 
colliding R-particles become infinite. However, there exist some possibilities to avoid the 
difficulty of superluminar velocities. These are discussed throughout the paper, together 
with other problems related to the correspondence of the elaborating scheme with the 
principles of Special Relativity. 

The organization of the paper is as follows. In Sect.2 we describe the principal 
system of the light cone equation and that of a “polynomial” worldline in “natural 
parametrization” Ao(t) = r and seek for the properties of its roots-particlcs. 
Remarkably, the (conserved) total energy (rest mass) of the system of R-C particles can 
take only integer values (equal to the higher degree of the parameterizing polynomials). 
One can also find therein a simple explicit formula that relates the value of the 
(conserved) total angular momentum to the leading coefficients of the parameterizing 
polynomials. 

In Sect.3 we investigate the other, so-called “timelike” class of polynomial 
worldlines, with p > n, where p and n are the higher degrees of the polynomials 
parameterizing the temporal and the spatial coordinates on the worldline, respectively. 
In doing this, we obtain the manifestly Lorentz invariant form of all three canonical 
conservation laws and discuss their physical meaning and interrelations with SR. The 
total energy (rest mass) appears again to be “self-quantized” being always equal to the 
number of roots. 

In Sect.4 we consider the (unexpectedly nontrivial) asymptotic behavior of the 
ensemble of R-C particles, at large values of the observer’s time T and for the class of 
timelike worldlines (defined by polynomials of a rather great order p). We demonstrate 
that, on the background of general runaway of particles, there are at least two time 
scales (Ti and To 3> T\) at which peculiar phenomena of “self-organization” of the 
particles’ ensemble do take place. Specifically, at T ~ T x separate R-C particles start to 
couple while at T ~ T 0 the pairs form themselves “multi-particle” clusters. Particular 
details of the process, including its time (ir)-reversibility, depend on the parity of p and 
n and their multiplicity. Finally, we give a graphical illustration of different stages of the 
process and speculate on the possible physical interpretation of the particles’ evolution 
(in particular, as a model of (in)clastic scattering). 

Concluding remarks concerning the advantages, drawbacks and perspective 
investigations of the “polynomial mechanics” are presented in Sect.5. In the Appendix 
A we prove (analyzing the structure of the so-called Sylvester matrices) explicit formulas 
for the total rest energy and angular momentum of an arbitrary system of R-C particles. 
In the Appendix B we present an illustrative example of collective dynamics of the 48 
(R-C) particles on a single polynomial worldline, calculate the exact values of all of the 
related conserved quantities and trace the most important stages of the evolution. 
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2. Conservation laws for polynomial worldlines in natural parametrization 

Consider a parametrically defined worldline X fl = F /t (r), fi = 0,1, 2, 3, with r being a 
monotonically increasing timelike parameter (not necessarily the proper time one!). The 
worldline itself defines the motion of one and only one pointlike particle. Nonetheless, 
if we consider the process of “detection” of the latter by an “observer” moving along 
its own worldline = /^(T), we encounter the effect of “multiplication” of the initial 
particle and formation of apparent copies of it located at different points on one and the 
same worldline. For this, one should make use of the well-known equation of the light 
cone of the observer (“the retardation equation”) [^j 

3 

(A'„(r) - &(T)) 2 = ~ £,(T)) 2 . (2) 

a= 1 

Generally, for any fixed position of the observer specified by the value of T, equation 
(J2| could have a lot of roots {Tj}, including real-valued ones, which define, with respect 
to the observer, the positions and collective dynamics of particles that all belong to the 
same worldline. However, it is well-known (see, e.g., [10]) that for the class of worldlines 
which contain no segments with superluminar velocities this is impossible: equation 
([2| has in this case only one root for which X 0 < £ 0 (the retarded solution) and one 
“unphysical” with X 0 > £ 0 (the advanced solution). 

Effect of “multiplication” for the class of purely tachionic worldlines and some of 
its interesting consequences were considered in the papers of B.M. Bolotovskii, V.L. 
Ginzburg and V.P. Bykov 0 0- as well as in the works of M. Ibison (see, e.g., 0). 
For the case of complex-valued worldlines an analogous construction had been presented 
in our works M. Remarkable properties of the worldlines on the complexified space- 
time background concerning, in particular, the procedure of generation of the Kerr-type 
solutions to the Einstein-Maxwell electro-vacuum system of equations, had been also 
numerously investigated by E.T. Newman mm, A.Ya. Burinskii mm et al. 

Generally, the class of worldlines containing segments with superluminar velocities 
and forbidden, therefore, in the canonical Special Relativity had been first proposed 
by E.C.G. Stueckelberg [I5j [16] who used these, in particular, to describe the 
annihilation/creation processes. Later such worldlines were examined by J.A. Wheeler 
and R. Feynman m h in the framework of the so-called “one-electron-Universe” 
paradigm. Of course, many problems arising from the permit of superluminar velocities 
are still on the agenda. However, there are reassuring expectations to avoid these 
difficulties: one of these will be considered in Sect.4. 

Our main goal here and below is to demonstrate the following remarkable fact. 
For any worldline X i M (r) for which the parameterizing functions F^{r) are arbitrary 
polynomials, corresponding set of roots-particles fixed by an “inertial” observer, obeys 
a number of conservation laws. The latter are defined for a set of SO( 3, l)-covariant 
quantities related to the roots of ([2]) and their time derivatives. The induced conservative 
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Note that we do not assume any Minkowski structure a priori; velocity of light is set to be unit, c = 1 
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dynamics is Lorentz invariant though different from the canonical relativistic one |^j 
Consider now the defining equations of an arbitrary “polynomial” worldline X /( (t) 
together with the light cone equation for an inertially moving observer + v y T. 

For simplicity, let us first accept the natural parametrization of the worldline X$ (r) = r 
and, as for the observer, assume him to be resting at the origin, £ 0 = T, £ a = 0. Then 
(J2| takes the form 

F := X(r) 2 + Y(t) 2 + Z(r) 2 - (T - r) 2 = 0, (3) 

in which 

[ X 1 = X(t) = a x r n + b x r n ~ l + ... + e x , 

< X -2 = Y(t) = a y r n + b y T n ~ l + ... + e y , (4) 

[ X 3 = Z(t) = a z r n + b z T n ~ x + ... + e g , 

are arbitrary polynomials of a highest degree n in r (some of the three polynomials may 
be of smaller degrees). Below we assume n > 2 and note that the leading term of degree 
2 n in (|3]) does not ever cancel. 

Then, at any instant of the observer’s time T, @> has exactly 2 n roots some of 
which are real while others enter in complex conjugate pairs. Corresponding copies- 
particles will be called R- or C-ones, respectively. The latter, as it was explained in the 
introduction, are located off the real worldline points but can be nevertheless visualized 
according to equal real parts of complex conjugate roots. 

Specifically, polynomial equation (J3| acquires the following structure: 

F = Ar 2n + Br 2n ~ x + Cr 2n ~ 2 + ... + (£> + 2 T)t + (E - T 2 ) = 0, (5) 

where the coefficients A, B, C ,..., D, E do not depend on the time T at all; in particular, 

A = a 2 +a 2 + a 2 = |a| 2 , B = 2 (a x b x + a y b y + a z b z ) = 2a ■ 6,..., , , 

£ = e 2 +e 2 + e 2 = |e1 2 . 1 j 

Now, making use of the first, linear in roots Vieta formula applied to the polynomial 
equation ([5]), we immediately obtain 

= —B/A = constant , (7) 

where summation here and below runs over all the roots i = 1,2,..., 2 n. Thus, we 
obtain the first “conservation law” for the quantity which is usually associated with the 
so-called “retarded time”. However, the values {r*} of the particles’ timclike parameter r 
may here be complex, decrease in the course of the monotonic increasing of the observer’s 
time T (then the corresponding root represents an antiparticle, see Sect.4), etc. In fact, 
only the latter “macro-time” is a well-defined evolution parameter of the theory. 
Consequently, the derivative quantity f := dr/clT will be conserved as well, 

5 > = 0 ' ( 8 ) 


3 For example, the field degrees of freedom do not contribute explicitly to the conservation laws; 
we suspect that these are indirectly represented in the characteristics of particles themselves, see the 
discussion below 
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while the next, quadratic in roots Vieta-like formula [^] arid its derivatives w.r.t. T lead 
to the following relations: 

t 2 = ( B/A ) 2 — 2 C/A = constant, TjTj = 0, r 2 + fjTj = 0.(9) 

Consider now one of polynomial equations for spatial coordinates of R-C particles, 
say for x = Xp 

H x :— x X (t) = x:— (a x T n + b x T n ~ l + ... + e x ) = 0. (10) 

Together with (|3]) it forms a joint system for determination of the implicit dependence 
of x on T ; to do this, one should eliminate r by taking corresponding resultant 
D x := Res[F, H x , r] and equating it to zero. Analyzing the structure of the determinant 
D x of the arising Sylvester matrix (see Appendix A for details), one obtains that the 
leading terms in the resultant equation will have the following structure: 

D x = Px 2n + Q x x 2n ~ l + (R x + E X T + G x T 2 )x 2n ~ 2 + ... = 0, (11) 


where P, Q x , R x , E x , G x are time-independent constants; in particular, we are interested 
below in the following two coefficients: 


P = A n =\a 


2 n 


G x = —nazA n = —naila 


\2n—2 


( 12 ) 


in which the quantity A was defined in (J6|. Quite analogous expressions for two other 
coordinates one will obviously come to. 

Now, making use again of the first two Vieta-like formulas, one obtains the following 
constraints between spatial coordinates of particles: 


E Qx , 

Xi = —— = const , 


P 


E 


U = (%) 2 - 2: 


.Rx + E X T + G X T 2 


(13) 


(14) 


P ' P 

and similar formulas for y- and 0 - coordinates of particles. From the linear constraints 


like (13) we see that the selected reference frame of the observer is just the center-of- 


mass one (on account of the assumed identity of R-particles which, therefore, possess 
equal masses) so that the total momentum is conserved and equal to zero'. 

= ^2(x' a )iP = 0, a = 1,2,3 (15) 

where / denotes differentiation w.r.t. r. 

Remarkably, conservation laws ( 13fT4 ) (as well as others below examined) are valid 
only for an inertially moving observer (including the currently considered case of the 
observer at rest). For nonlinear polynomial worldlines of an observer (corresponding 
to a non-inertial reference frame) the structure of resultants like ( [IT] ) becomes more 
complicated, and conservation laws turn out to be broken. 


4 Precisely, we make use of the formula for the sum of squares of the roots which can be obtained from 
the quadratic Vieta formula with account of the first, linear in roots one 
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Composing now the SO (3)-invariant combination from the sums of squares like (14) 
and differentiating twice the result by T, one obtains the following conservation law for 
the ensemble of R-C particles: 


W := ^2(x a )i(x a )i + ^2{x a )i(x a )i = -2 


G x + Gy + G z 
P 


= constant. (16) 


Substituting then expressions ( |12[ ) for the constants P, G x (and for corresponding ones 
defining G y , G z ) we obtain finally 


W = 2 n- 


ai + a,. + a; 


A 


= 2 n. 


(17) 


Thus, we have obtained a sort of quantization rule for the quantity which, in fact, is the 
analogue of total energy of the R-C ensemble: 


W ir 

— = K + V = const = n, 


(18) 


where the first term K := | J2(.%a)i(x a )i reproduces the total kinetic energy while the 
second one V := \ J2(^a)i(x a )i stands for the quantity known in classical mechanics as 
virial and proportional (for radial forces homogeneously depending on mutual distances) 
to the total potential energy (see, e.g., m FT8]). Consequently, any “polynomial 
worldline” (in natural parametrization X 0 (r) = t) of a highest degree n defines a 
conserved energy-like SO (3) -scalar quantity which is moreover universal, integer and 
equal to the degree n. 

Let us demonstrate further that the vector of total angular momentum 

M = f x r (19) 

is also conserved for any R-C system of particles. To do this, implicitly differentiating 
the principal equation (J3|, we come to the rational expression for t (in which F : = 
dF/dT, F':=dF/dr): 

F T -r 

r = ~¥' = l(fGGf)^rxx r TYY r +Z~Z' ^ 

and then obtain for a component of angular momentum of a particle, say for M x , 
expression of the form 

M x = YZ - ZY = ( YZ' - ZY')t. (21) 


On account of (20), expression (21) reduces (except in the singular points of mergings) 
to the following equation: 

M X F' + (YZ 1 — ZY')F = 0 ^ 

N x := M X {{T - t) + XX' + YY' + ZZ') - (T - t)(YZ' - ZY') = 0, 1 } 


in which the polynomial N x (t , T, M x ) is of the degree 2 n — 1 in r (and linear in T and 
M x ). Eliminating then r via taking the resultant R x := Res[F, N x ,r] and equating it 
to zero, we come to a polynomial equation R X (T,M X ) = 0 of the expected degree 2 n in 
M x . The latter implicitly defines the dependence of (the first component of) angular 
momentum of particles on the observer’s time T. 








By virtue of the specific structure of the determinant of the corresponding Sylvester 
matrix (see the Appendix A), we obtain that the two leading coefficients in the resultant 
equation are of the following appropriate form: 

R x = A (T)Ml n + A {T)a x Ml n ~ l + ... = 0, (23) 


where a x is a time-independent constant (see the Appendix A) while A (T) is a specific 
polynomial in T, the discriminant of the principal equation ([3]) , which turns to zero 


only at singular points of mergings. The above obtained structure of resultant (23), 


together with the first, linear in roots Vieta formula, guarantees the conservation of 
total quantity M x : 

^( M x )i = -a x = const. (24) 

i 

The same, of course, is valid for two other components of the vector of total angular 
momentum M. 

To corroborate the property of conservation, we derive in the Appendix A the 
following universal explicit formula for the value of M: 

a x b 


M = 2- 


(25) 


which turns out to depend only on the leading coefficients a = {a x , a y ,a z }, b = 
{b x ,b y ,b z } of respective degrees n and n — 1 in the defining polynomials X a (r), a = 
1,2,3. If only one polynomial, say X\ = X f(r), is of a degree n while the two others 
have lower degrees, the total angular momentum is identically zero. 

Finally, we announce that for the more general parametrization Xo(r) = S{r) ^ r 
for which the degree p of the polynomial S'(r) is smaller than n, that is 0 < p < n 
(the so-called “spacelike” worldlines, see below), the above conservation laws all remain 
valid. However, the proof of this statement as well as corresponding explicit formulas for 
the total angular momentum (25) and the analogue of total energy © become much 
more complicated and depend on the particular degrees of generating polynomials Q. 


3. “Timelike” polynomial worldlines and Lorentz invariant dynamics 

The principal retardation equation (|3]) , together with the defining equations of a 
polynomial worldline Q, compose by themselves a Lorentz-invariant system. On the 
other hand, we recall that only for an inertial worldline of the observer £ m (t) the 
conservation laws for the system of R-C particles remain valid. In that way, we encounter 
a remarkable relation between the (quadratic) structure of the light cone equation and 
the conservative character of the collective dynamics it induces in an inertial reference 
frame. 

Specifically, under a boost, the parameters r and T in (J3]) should be considered 
invariant while the four coordinates of both the observer £ /t (T) and particles X y (r) 
transform in the canonical way. This means that in a new inertial reference frame the 
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above used simplest gauge A 0 (r) = r will be broken, and the explicit formulas for 


conservative quantities like (17) or (25) will be no longer valid 


Thus, only the SO(3, l)-covariant combinations of conservative quantities should 
be considered as physical meaningful. Such combinations, on the other hand, can be 
constructed only from the constraints related to the linear or quadratic in roots Vieta- 
like formulas [^j As for the first one, the corresponding quantities A ;i = {Ao, X a } 
evidently form a 4-vector. Since its spatial part, on account of equal (unit) masses of 
the roots-particles, had been identified with the momentum , the 4-vector A M as a whole 
should be treated as the energy-momentum 4-vector of an individual particle. 

However, for the “observer at rest” in the chosen gauge A 0 (r) = r from Eqs. 


and (15) it immediately follows that the total energy as well as the total momentum 


of the whole R-C ensemble is identically null; this, of course, will then be true in any 
inertial reference frame and cause serious difficulties concerning the identification of R- 
C particles in the case of the above used parametrization A 0 (t) = r and, generally, in 
arbitrary case of the “spacelike” worldlines (i.e., for A 0 = S(t), with deg(S) = n < p ). 
Fortunately, the situation changes drastically for the class of “timelike” polynomial 
worldlines which we are going to examine below. 

Specifically, let the parameterizing polynomials of the worldline under consideration 
have the form (cf. Q): 


X 0 = S(t ) = a s r p + b s r p 1 + ... + e s , 
X\ = X(t) = a x r n + b x r n ~ l + ... + e x , 

^n— 1 


(26) 


X 2 — Y (t) = a y r n + b y T n 1 + ... + e y , 

X 3 = Z(t) = a z r n + b z T n ~ l + ... + e z , 

where we assume p > n > 2 and a s ^ 0. Then, for sufficiently large values of r, 

corresponding values of the (proper) Minkowski intervals of R-particles ds 2 = gLAq — 

dXl — dX^ — dXl ~ p 2 a 2 T 2p ~ 2 dT 2 will be evidently positive definite ; this justifies the 
terms “timelike” or, vice versa, “spacclike” (precisely, asymptotically timelike/spacelike) 
for the characterization of the two different classes of polynomial worldlines. Needless 
to say that for timelike worldlines the natural parametrization Aq(t) = r is forbidden; 
instead, the “polynomial reparametrization” r = Ajfi + Bj-ifX 1 + ... is allowed that 
preserves the property of the initial worldline to be (asymptotically) timelike. 

For a timelike worldline (26), the principal equation of the light cone (j2]), for 
observer at rest, acquires the form (cf. 


an 


F \= (T — S(t )) 2 - A(r) 2 - Y{t) 2 - Z{t ) 2 = 0, 


(27) 


so that the polynomial F is of the degree 2 p in r, and at any T the observer detects 2 p 
associated R-C particles, ft is easy to see that all the conservation constraints like ([T] 
[9]), for the roots r*, % — 1,2,..., 2p or their time derivatives hold equally in the timelike 
case. As regards for the corresponding timelike coordinates s* of individual particles, 

5 E.g., formulas for the sums of cubes of the roots or their time derivatives do not possess well-defined 
tensor properties 
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one can make use of the auxiliary equation 
G s — S(t) = 0 


to simplify the equation of the light cone (27) to the following form: 
F = Q 2 -X 2 - Y 2 -Z 2 = 0, 


(28) 


(29) 


in which the quantity Q := s — T does not contain the parameter r so that (29) is now 


of a smaller degree n in r. Now the structure of resultant of two polynomials (28) and 


(29), that is of R s := Res[G, F,t], becomes quite transparent, and after eliminating r 
one obtains the resultant equation R s = 0 in the following form: 

R s = Q 2p + a{s)Q 2p ~ 2 + (3(s)Q 2p ~ 4 + ... = 0, (Q = s-T ), (30) 

where a, (3,... are some polynomials in s of degrees at most 1,2,..., respectively. 


Collecting then in (30) the leading terms in s of highest degrees in T which are, in 


fact, contained only in the first term Q 2p , one obtains 


R a = s 2p - 2 pTs 


2p— 1 


+ 


2p{2p — 1) 


T 2 s 2p ~ 2 + ... = 0 


(31) 


(note that all other terms do not enter into the principal characteristics written out 
below and can be discarded). 

Making now use of the corresponding Vieta-like formulas and taking their time 


derivatives, one obtains from (31), in full analogy with the previously described 
procedure (Sect.2): 

E := ^2 Si = 2p = constant, (32) 

W 0 := + X/ $ iSi = - 2—— = 2 p = const. (33) 

As to the total momentum P = {P a }, a = 1, 2, 3 and conjugate conserved quantities 
{W a } related to the quadratic in roots Vieta-like formulas, following the same procedure 
one immediately concludes that these both are null in the timclike case, 

Pa := X^a)i = 0, (34) 

and 


W a := y 2(ig) 2 i + = 0. 


(35) 


We see therefore that, in the timelike case, identification of the quantity P t , := X fJ = 
{X 0 , X a } with the energy-momentum 4-vector becomes quite adequate. Indeed, the 
observer at rest (when the reference frame coincides with the center- of-mass frame of 
the R-C particles’ ensemble) will get the positive-definite total energy E which, according 


to (32), is moreover integer, universal and equal to the number of roots-particles 2p! One 


can thus introduce the SO (3, l)-invariant 

M 2 := E 2 -P 2 = (J^Si) 2 - 

which should be evidently identified with the square of total rest mass of the particles’ 
ensemble. 




— 4p 2 > 0, 


(36) 
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Quite analogously, from the conserved quantities (33),(35) one can construct the 
SO(3, l)-scalar 

W := (W 0 — W 1 - W 2 - W 3 ) = 2p = M > 0 (37) 

which turns out to be precisely equal to the positive-definite total rest mass of the R-C 
system! Note that in the 2D non-relativistic case considered in corresponding 

SO (3)-invariant quantity conserved by virtue of the quadratic Vieta-like formulas also 
represents the analogue of total mechanical energy of the system; in the Lorentz invariant 
dynamics under consideration both “definitions” of the total rest energy (rest mass) do 
coincide, and we regard this fact as indicative of the self-consistency of the theory. 

On the other hand, the structure of quadratic relations like (35) demonstrates that 
the ordinary kinetic-like part of the total energy represented by the first term therein 
would be constant only under account of the second term, the (relativistic analogue of) 
potential-like part. The latter, however, directly depends on accelerations of particles 
and, in a sense, substitutes for the energy of the fields associated with particles in the 
canonical approach. We shall return to discuss these issues in the conclusion. 

It is also noteworthy in this connection that in the rest frame the equality of two 
quantities (32) and (33) follows, in fact, from the following identity for the roots of an 
arbitrary timelike system of equations (26),(27) and their time derivatives 

= + 




(38) 


Deep physical sense (if any) and mathematical origins of this remarkable identity related 
evidently to the total rest energy (rest mass) of the R-C system is not yet revealed. 

To conclude, let us write out the manifestly Lorentz invariant form of the 
conservation laws arising in the scheme under consideration. For the energy-momentum 
4-vector one gets 




dT 


— Ca — const, 


while for the scalar of total mass one obtains 


E 


dX,#) dX (t) 

dT dT 


d 2 Xf ) 

E X / J «^W = const 


(39) 


(40) 


Note that, instead of the individual proper times of the constituent particles in the 
canonical approach, one deals here with the unique Lorentz invariant proper time of the 
(inertially moving) observer T. 

Let us finally clarify the situation with the total angular momentum. Recall that in 
a relativistic invariant scheme one should consider the skew symmetric tensor of angular 
momentum with six independent components 

M \^i = E XfoUv], (41) 

where the quantities LQ are the 4-velocities of particles. However, for the considered 
R-C system, the latter should be defined via differentiations w.r.t. the observer’s proper 
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time T, that is 


TT - d ^ = X 

U>1 ~ dT “ A " 


(42) 


Three spatial components M a := ^£ a bc.M[bc] of the tensor (41) evidently reproduce the 
ordinary expression for the angular momentum vector 

M a = £ a bcXbX c , (43) 

while the components M[ 0a ] are their relativistic counterparts. In our scheme, all these 
six quantities are conserved for the system of R-C particles as a whole though, in the 
timelike case, we were unable to find explicit formulas for them resembling (25). The 
SO(3, 1)-invariant corresponding to the conserved tensor (41) is evidently its 4-square 
Contrary to the square of the energy-momentum 4-vector P y = X y the 
latter does not, however, take a universal value. 


4. Asymptotic behavior: formation of clusters &; scattering pattern 


Collective dynamics of R-C particles, in the timelike case, reveals universal remarkable 
properties at large values of the observer’s time T. In fact, asymptotically, individual 
particles join into compact groups, clusters, which contain two or four R-particles and an 
essential number of C-particles each. At the same time, different clusters move off from 
the origin with decreasing velocities. Specific details depend on the particular values of 


the degrees p and n of the generating polynomials (26), and we briefly examine them 
below. 

At sufficiently large values of the observer’s time T, the terms of highest degree 


2 p dominate in the light cone equation (27) parametrized by (26). Therefore, equation 


(27) asymptotically takes the obvious form: 
(T - a s r p f ~ 0, 


(44) 


with the following, twice degenerate set of solutions (to be concrete, let us first assume 
a s > 0 and T > 0): 

r k ~ |T/a s | 1/p exp (2mk/p), k = 1, 2,... ,p ,..., 2p. (45) 

For the limiting values of the (real parts of) coordinates of particles we obtain then from 


(26): 


{X c }k — 3?{a c r^} = a c \T / a s \ n ^ p cos (2nkn/p), c— 1,2,3. (46) 

(here and below we rename {a x ,a y ,a z } = {a c } = { 01 , 02 , 03 } and so on for { 6 C },...). 
However, if for some values of the angle := 2nk(n/p) the coordinates {X c }k turn 
to zero, one should take into account the next terms of the degree n — 1 , so that the 
particles are located much closer to the origin and move slower, respectively. We can 
call such set of particles exceptional. 


The asymptotic form (44) of the principal light cone equation and corresponding 
expressions (46) for coordinates of particles hold for T T 0 , T 0 = a s {h s /a s ) p . The 
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first order corrections due to the account of the terms of the degree 2p — 1 lead, instead 
of (44), to the following approximate form of the light cone equation: 

(T - a s r p - b s r p ~ l f ~ 0, (47) 


valid at the smaller time scale T 0 > T ~ Ti, where Tj = a s {\a\/a s ) p ^ p ~ n \ Finally, the 
second order corrections take into account the highest order terms (of degree n ) of the 
polynomials {A c } and lead to the following representation of the principal equation: 

(T — a s r p — b s T p ~ 1 ) 2 — \a\ 2 r 2n ~ 0, |a | 2 := a\ + al + a 2 , (48) 


which is valid at T ~ Tj and breaks the degenerate structure of the roots-particles 
related to the full square in the right hand side of (47) or (44). 

ft seems quite natural (and most interesting from the physical viewpoint) to 
examine the case when the whole number of roots-particles N = 2p is great, p > 1. 
As for n, it defines, as a rule, the number of clusters and the number of particles in a 
cluster which, evidently, are in an inverse proportion to each other. Elementary analysis 
of the successive approximations (48), (47) and (44) of the exact light cone equation 
(27) leads then to the following conclusions. 


A. At the interval 0 < T < Tj . 

1. At particular instants of T, some two of R-particles merge (annihilate) and 
transform themselves into a composite C-particle or vice versa. Before an annihilation 
event, the course of “retarded” time r reproduces that of the observer’s time T for one 
of the merging particles and is inverse for the other one. Thus, one can regard the 
first/second formation as a particle/antiparticle, respectively. 

2. At T ~ T\ all the “events” are over and, in what follows, the whole number 
of R- and C-particles remains invariable. Specifically, there are finally only two (in the 
case of odd p) or four (for even p) R-particles. The number of C-particles may be thus 
very large. This early stage of evolution for the case p = 24, n = 20, (A = 2 p = 48) is 
depicted at figure [l] 


B. At the interval T\ -C T < T 0 . 

3. All R-particles break up into one (for odd p) or two (for even p) contracting 
pairs (R-pairs) while all C-particles form C-pairs consisting of two slightly separated 
and approaching each other pairs of complex conjugate roots (that is, of four roots- 
particles altogether). Thus, one observes the phenomenon which can be called the “first 
phase transition" (see the illustration for the case p = 24, n = 20 at figure [2]). 

C. At the interval T T 0 . 

4. All (non-exceptional) particles-roots concentrate in vicinity of the straight line 
defined by the ratio of the coefficients a x : a y : a z and undergo a runaway according to 
the law cx T n / p , with decreasing velocities and negative accelerations (note, however, 
that for p n the motion is nearly uniform). As to exceptional C-pairs, they move 
along the other direction specified by the ratio of coefficients b x : b y : b z \ their recession 
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Figure 1 . First stage of evolution, T = 10, T <C T\\ R- (in red crosshairs) and 44 
C-particles (in green points), some of these being off the depiction space; the pairs are 
not yet formed 


law is oc T( n ~ 1 " p so that they are located much closer to the origin and move slower, 
respectively. 

5. For multiple p and n, formation of clusters takes place (the “second phase 
transition ”); one of these consists of either one (for odd p) or two (for even p) R-pairs 
“dressed” by a number of closely disposed C-pairs. All other clusters consist solely of 
C-pairs. For p = 24, n = 20 this last stage of evolution is depicted at figure [3]|^j The 
whole number of clusters, of pairs (as well as particular configurations of the latter 
within a cluster) depend on the parity and multiplicity of p and n. In the course of time, 
separation of clusters evidently grows whereas mutual distances between pairs within a 
cluster (as well as between particles constituting any pair) relatively (i.e., in comparison 
with distances between clusters) decrease. 

D. “Full history” at the whole interval —oo > T > +oo . 

6 In the considered example there are no exceptional particles at all and, thus, only one spatial direction 
of runaway well noticeable at figure |3l 
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Figure 2. Second stage of evolution, T = 10 11 , To > T T\\ one R-pair (the other 
one is far off the frame) and 22 C-pairs; all the pairs are well distinguishable and start 
to gather into clusters. 


Making use of the same considerations, one can analogously describe the R-C 
dynamics on the negative semi-interval and relate, in particular, the incoming and 
outgoing distributions of particles. As a result, one comes to the following conclusions. 


6 . For odd p, the asymptotic equations (44) and (46) are invariant under the 
substitution T (->■ — T, r H» — r, X c H > ±X C , for even/odd n, respectively. Therefore, 
one has an “asymptotically time-symmetric” process resembling the picture of elastic 
scattering of a beam of clusters in the transmission/reflection regimes, respectively. 

7. For even p, equation (44) possesses the symmetry T i->- — T, r i—> ±?,r, X c 

— X c (for even n only). In particular, for —oo < T •C —To it has no real roots (so that 
only C-clusters are present in the incoming beam) while at +oo > T T 0 it possesses 
4 real roots corresponding to a single R-cluster (either “dressed” or “naked”, see also 
item 5). Thus, the process is “asymptotically time-asymmetric” and model that of an 
inelastic scattering. At small values of time T one encounters a cascade of creations of 
R-pairs which successively annihilate then into composite C-particles; the final stage of 
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Figure 3. Third stage of evolution, T = 10 50 , T Tq: four well-formed clusters run 
away along a distinguished 3-direction; two R-pairs are designated by red circles and 
22 C-pairs by green squares; individual particles constituting pairs are too close to be 
seen. 


the process (with one “dressed” R-cluster and a lot of C-clusters) was described above, 
see again the illustration at figure |3j 

5. Conclusion 

In fact, we have demonstrated that a quite nontrivial and well defined collective 
dynamics of pointlike particles can be induced solely by generating functions defining 
a single worldline. For all this, we did not make use of Lagrange functions, systems 
of differential equations or any other ingredients of canonical mechanics. When the 
worldline is defined in an ordinary parametric form, an apparent set of particles arises 
as the set of roots of the equation of light cone of an observer moving along its own 
worldline. 

We investigated the simplest case of a “polynomially parametrized” worldline. Then 
it is necessary to consider two different kinds of copies-particles, R- or C-ones, that 
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correspond to real or complex conjugate roots of the light cone equation, respectively. 
Then all physical characteristics should be taken as functions of the proper time T of 
the observer. 

Remarkably, in the reference frame of an inertially moving observer (and only 
in that) the induced collective dynamics turns out to be conservative , that is, for 
arbitrary polynomial worldline obeys a number of conservation laws. Their structure is 
completely governed by the Vieta formulas for the roots-particles and does not include 
thus any characteristics of the fields produced by the particles but only the particles’ 
characteristics themselves. On the other hand, the presented scheme is based on the 
effect of retardation and is manifestly Lorentz invariant. Moreover, the energy-like 
dynamical characteristics contain the term that depends on the accelerations of particles 
and can, in principle, substitute for the held degrees of freedom (see below). 

As to the (electromagnetic, curvature, etc.) fields themselves, they are uniquely 
defined via the principal spinor of a shear-free null congruence of lightlike rays produced 
by pointlike objects moving along an arbitrary single worldline on the real space-time 
or even on its complex extension (see, e.g., munununii). Moreover, at the instants 
of merging of a pair of particles the fields form a singular null straight line, caustic , 
that connects the points of merging and observation and is, perhaps, responsible for the 
information transmission P|. However, the energy and momentum of the fields, as we 
have seen, do not explicitly contribute into the structure of conservation laws. 

We have shown that there exist two essentially different classes of polynomial 
worldlines, the so-called “spacclike” or “timelike” ones Q and the latter seems to be 
more suitable for physical interpretation. For the whole R-C system of particles (or 
an approximately isolated subsystem of these), in the timelike case, all the formulas of 
the canonical relativistic mechanics are valid, including the mass shell equation (B.4). 
However, expressions for the momentum or energy of an individual particle preserve 
in fact their non-relativistic form and, taken separately, do not satisfy the mass shell 
constraint. Moreover, velocities of individual particles can exceed the light limit and 
infinitely increase before their merging (annihilation of a pair). However, at late stages 
of the evolution, when all the mergings (“events”) are over, velocities of all the particles 
decrease and are certainly less than that of the light (Sect.4). 

Anyway, the algebraic dynamics under consideration, though relativistic invariant, 
strongly differs from the canonical mechanics of Special Relativity. On the other hand, 
it is well known that there are a lot of “eternal” problems in the application of the 
SR-formalism itself (contrary to that of the non-relativistic mechanics). For example, 
the energy of fields diverge for the pointlike objects; physical characteristics of particles 
themselves are, in fact, compared at different proper time instants ( different retardation 
intervals), etc. In this connection, there were a number of attempts (see, e.g., jTSJjZO, 21]) 
to reformulate the theory in the terms of the unique (laboratory) time. To do this, 
the authors proved it possible to totally substitute the field degrees of freedom by the 


7 The third, rather remarkable “lightlike” class of worldlines corresponding to the case p = n is under 
study 
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infinite set of higher order time derivatives of the coordinates of individual particles. 
The situation arising in the above presented scheme seems to be rather close to this 
picture. 

As to the most interesting properties of the developed algebraic dynamics, these 
are, of course, the “self-quantization” of total rest energy of the R-C system of particles 
(which can take only integer values) and the general recession of particles’ ensemble. 
The latter is moreover accompanied by formation of pairs and then clusters of particles 
at very large values of the observer’s proper time T T 0 . Qualitatively, the arising 
picture resembles that of (elastic/inelastic) scattering of a beam of (composite) particles. 
One can even speculate about some likeness between the scenario of evolution of the 
system of clusters and the cosmological evolution of the Universe as a whole. 

As an obvious generalization of the above presented construction, one can consider 
the class of worldlines parametrized by rational functions. Preliminary investigations 
demonstrate that for a wide class of rational worldlines the whole set of conservation 
laws will be satisfied as before. However, the asymptotic behavior in this case will be 
quite different and even more remarkable from the physical viewpoint. 

Among the obvious physical drawbacks of the theory one can note the effect of ‘de- 
isotropization of matter’ which takes place at large values of the “cosmological” time 
T: all the particles run away along one or two privileged spatial directions (see item 
4 in Sec. 4). Unfortunately, there is also a negligibly small number of R-particles, at 
least asymptotically. However, the situation could become much more realistic in the 
case of an implicitly defined polynomial worldline [310! which, in a sense, substitutes 
for the whole set of ordinary, parametrically defined ones. Corresponding construction 
(exploiting again the light cone equation) will be considered elsewhere. 

In any case, above revealed purely mathematical properties of the fundamental 
equation of the (polynomially parametrized) light cone are quite unexpected and seem 
to be remarkable by themselves. On the other hand, they explicitly relate to the well- 
known fundamental properties of physical matter (conservation, quantization, formation 
of clusters, recession, etc.). We believe that this correspondence is not casual, and a 
lot of other, physically relevant, “numerical” relations are encoded in the mathematical 
structures like those examined in the paper. 


Appendix A. Structure of Sylvester matrices and conservation laws 


Let us proof the structure © of the resultant D x = Res[F, H x: r] of two polynomials 
F{t) and F[ x = x — X(r). We start from the simplest case n = 2; then the two 
corresponding equations (|5]) and ( 10 ) are of the form 


F = At 4 + Bt 3 + Ct 2 + (D + 2T)r + (E - T 2 ) = 0, (A = \a\ 2 ± 0), 
and 

H x = a x r 2 + h x r + (e x - x) = 0 , (a x 7 ^ 0 ). 


(A.l) 

(A. 2 ) 
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The 6x6 Sylvester matrix E x for the above resultant will have the following form: 


Eh = 


( A 
0 
0 
0 
0 

\ &X 


B 

A 

0 

0 

K 


C 

B 

0 

&X 

bx 

e x — x 


D + 2T 
C 

b x 

e a 

0 


x 


E-T 2 
D + 2T 
b x 

0 
0 


x 


0 \ 

E-T 2 

e x — x 

0 
0 
0 


(A.3) 


To obtain the resultant, one should compute the determinant D x of (A.3) which is 
obviously a polynomial of the degree 2n = 4 in x; we, however, are interested mostly 
in three leading terms of degrees 4,3,2 in x written out in ( fTTj ). Their structure and, 
in particular, (in)dependence on the time T follow directly from the form of Sylvester 
matrix (A.3) and completely correspond to that adduced in the text: 

D x = det E x = Px A + Q x x 3 + (R x + E X T + G x T 2 )x 2 + ..., (A.4) 

where, as it is easy to see, P = A 2 ,...,G X = —2 Aa 2 . Thus, the resultant equation 
D x = 0 leads to the conservation laws of the form represented by (15) and (16). It 
is evident (and had been confirmed by the computer algebra computations) that a 
quite analogous structure of the above determinant-resultant will be observed for any 
n. Specifically, for the leading coefficients P and G x one obtains the expressions (12) □ 
Let us examine now the structure of resultant R x := Res[F, N x ,r] of the two 
polynomials F(t) and N x (where the latter is represented by (22) and related to the 
x-component of the angular momentum (21)). For this case the Sylvester matrix has 
the form essentially different from (A.3). Specifically, for n = 2 the principal light cone 
equation is represented by (A.l) while the polynomial equation (22), on account of the 
defining equations of the worldline Q, takes the form 

N x = M x (4At 3 + 3 Bt 2 + 2 Ct + D + 2T) - (T - t)(S x t 2 + L x t + K x ) = 0, 
where, in particular, 

S x 2{bya z aybz') 
and all S X ,L X , K x are independent on T. 

Designating now for short D = D + 2T, E = E — T 2 and omitting for a time the 
component index x, we are able to write out the corresponding 7x7 Sylvester matrix 
of the resultant R = Res[F , N, r\. 

( A B 

0 A 

0 0 

0 0 

0 0 

0 4 AM + S 

\ 4AM + S 3 BM + U 

where S(— S x ) = const 
—TK x are linear functions of T. 


(A.5) 


(A.6) 


c 

D 

E 

0 

0 

\ 


B 

C 

D 

E 

0 



A 

B 

C 

D 

E 



0 

4 AM + S 

3 BM + U 

2 CM + V 

DM + W 

( 

A.7) 

4AM + S 

3BM + U 

2 CM + V 

DM + W 

0 



3 BM + U 

2 CM + V 

DM + W 

0 

0 



2CM + V 

DM + W 

0 

0 

0 

) 


while U(= U 

e) = L x — TS X , V(=V X 

5 

1 

W x ) = 
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Now it is easy to see that the first two leading terms in the determinant of (A.7) 


of degrees 4 and 3 in M do not involve any of time-dependent coefficients U, V, W. The 
leading term of degree 4 does not also contain S and, whether all these coefficients are 


set zero, the determinant of (A.7) becomes equal to M 4 Res[F , F', r] = M 4 A(T), that is, 
proportional to the discriminant A (T) of the principal light cone equation F — 0. The 
next term of degree 3 contains the constant S, is also proportional to the discriminant 
A(T) and equal exactly to M 3 A(T)(S/A) so that one has finally: 

R {n =i) = A(T)M 4 + aA(T)M 3 + ... = 0, ( a = S/A ). (A.8) 

It is quite easy to observe now that the generalization of the above considerations to any 


n > 2 leads to the structure of resultant R x = Res[F, N x ,r] similar to (A.8), namely 
(we can restore now the component index x): 

= S x /A , (A.9) 


R x = A (T)Mf + a x A{T)Ml n - 1 + ... = 0, 


a. 


that is, to the structure completely identical to that of (23) presented in the text. 
The quantity — a x which determines the (x-component of) the total vector of angular 
momentum is in fact equal to the ratio of the leading coefficients S x (= S ) and A of 


corresponding polynomials. Substituting their expressions (A.6) and (A.l) through the 
coefficients of defining polynomials Q one obtains finally: 

(a. io) 


al + a? y + a 2 ’ 


and, of course, analogous formulas for two other components. Thus, we proof the 


universal formula (25) for the vector M of total angular momentum of the system of 


R-C particles. It is noteworthy that the expression for M remains the same for any 
degree n of the defining polynomials Q. □ 


Appendix B. Dynamical characteristics for the examplified worldline 

In Sect.4 we presented a number of figures illustrating the collective dynamics of 2 p = 48 
roots-particles for the timelike (p > n) worldline with p = 24, n = 20. Specifically, we 


accepted the following (randomly selected) form of the polynomials (26) parameterizing 
the worldline: 

S = t 24 + 5r 23 — 2 t 22 + 21t 19 — r 7 + 5r — 8, 

X = t 20 — Hr 19 — 13r 18 — 3r + 5, 

Y = 2r 20 + 17 t 19 + 7r 18 + 2r 2 - 1, 

Z = —3t 20 + 13r 19 - 9 t 17 + 7r 4 + 12. 


(B.l) 


Substituting (B.l) into the light cone equation (27), taking the resultant Res[F,G,r 


of its left-hand-side polynomial with that of G represented by (28) and equating the 


resultant to zero, one obtains the following polynomial equation for determination the 
dependence s = s(T ): 


s 48 - (48T + 5)s 4/ + (1128T 2 + pT + ^)s 4b + ... = 0, 


(B-2) 
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where 6, ip , -0 are some great integer-valued real numbers (their exact values being 
irrelevant in what follows). From (B.2), for the sums of all 48 roots and their squares it 
follows immediately: 

J2 Si = (48T + 5), J2 S 1 = ( 48T + £) 2 - 2 • (1128T 2 + pT + if;), (B.3) 

i i 

so that after appropriate differentiations w.r.t. T we obtain the following conserved 
values for the rest energy (rest mass) of the system of R-C particles under consideration: 

E = ^2 Si = 48 = 2p, W 0 = ^2^ + SiSi = 48 = E, (B.4) 

i i 

in full agreement with the general formulas ( J32|33 ) or (38). Of course, it is easy to 
assure that the total momentum P represented by (34) as well as the corresponding 
quantities W a as given by (35) turn to zero in the assumed rest frame of the observer. 

Consider now the skew tensor of angular momentum (41). All its six components 
are nonzero and conserved but their direct calculation by means of the procedure 
described in Sect.2 and Appendix A is rather cumbersome since their values turn out 
to be unexpectedly great. Nonetheless, with the help of the computer algebra system 
Mathematica, we find in the end: 


M x = 5585145722802, M y = 11642674433710, M z = 316030077916, 

K x = 2277192433512854, K y = -799390472683014, K z = -4701753623402064, 


(B.5) 


where the components {M a } are defined by (43) while the quantities {K a } correspond 
to the components M[ 0a ] of the angular momentum tensor. The two 5'0(3)-invariants 
are, respectively, 

"|2 —. nj2 i nj 2 
>2 


f 2 + M 2 + M 2 

K 1 = Kl + K' 2 y + Kl = 27931117642239004403999706809608, 


M 2 := Ml + M 2 + M 2 Z = 16684559574445537237998360, 


(B.6) 


while the conserved Lorentz (SO(3, 1)) scalar is 

' K 2 — M 2 = ...?» 2.793111763 • 10 31 . 


(B.7) 


Of course, these exact results have been checked by direct numerical calculations of the 
sums of roots and respective dynamical quantities of these for different values of the 
observer’s time T. 

The discriminant polynomial equation A (T) = 0 possesses 8 real roots defining 
the instants T merge of the annihilation/creation events. Specifically, at the interval 
ti < T < t 2 , ti — —18.8923, t 2 — —18.3244 the light cone equation has 8 
(maximal possible number for the considered example) real and 20 pairs of complex 
conjugate roots defining 8 R-particles and 20 (composite) C-particles, respectively. At 
the moment T ~ t\ one has a creation of a pair of two R-particles (precisely, of a pair 
particle/antiparticle) while at the moment T ~ t 2 there is a merging of two R-particles 
accompanying by the formation of one C-particle, and so on. The first event (creation 
of the first R,-pair) takes place at the great negative value T ~ —0.7548374411 • 10 16 . 
The last event occurs at T ~ 4.434361164, and after this moment the number of R- and 
C-particles remain invariable (four R-particles and 22 C-particles). 
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The critical values T\ =~ 2.7 • 10 3 and T 0 ~ 6.0 • 10 16 as defined in the text (Sect.4) 
mark the intervals at which the particles are separated (0 < T < Tf), join into pairs 
(Ti < T < T 0 ) or start to gather into clusters (T T 0 ), see figs.2 and 3, respectively. 
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